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Abstract 
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ble curve of current supernova observations with a non-accelerating cosmological model. We also 
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to other speed of light communications such as gravitational radiation and neutrino fluxes. 
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I. INTRODUCTION 



The study of the intergalactic medium (IGM) has long been an important field in both 
astrophysics and cosmology. Current study of the influence of the IGM on distance red- 
shift relation is mainly focused on light absorption, e.g., the magnitudes of supernovae are 
corrected for light absorption by the IGM when drawing a Hubble diagram. We know that 
light paths in a dielectric medium are different from those in vacuum. According to classical 
electrodynamics light is altered in both speed and direction (light refraction). The impact 
of light refraction on the cosmological distance redshift relation is interesting from both a 
theoretical and observational point of view. An interesting and useful theoretical tool to 
study light refraction in curved spacetime is Gordon's optical metric The idea of the 
optical metric is simple, any solution to Maxwell's equations in a curved spacetime filled 
with a fluid whose electromagnetic properties can be described by a permittivity e(x) and a 
permeability fi(x) can be found by solving a slightly modified version of Maxwell's equations 
in a related spacetime with vacuum values for the permittivity and permeability, i.e., with 
e(x) = 1 and fi(x) = 1. 

Our motive for undertaking this work was two fold. First we wanted to revive the old 
optical metric theory of Gordon l|] and show how one applies it to a modern cosmological 
model. The Second was to join in current efforts to find alternative explanations of the 
recently observed cosmic acceleration besides the existence of an exotic p = —pc 2 material, 
i.e., besides the cosmological constant (see 

aaaasaaai Ifl,llll,ll2|,ll3|,ll4|). In Sec- 



tion [TT] we introduce the reader to Gordon's optical metric. In Section IHII we use the WKB 
approximation to derive the transport equations for the amplitude and polarization vector 
of waves moving through a spacetime possessing an optical metric. Readers not interested 
in the mathematical details can skip this section. In Section IIVI we construct the optical 
metric for Friedmann-Lemaitre-Robertson- Walker (FLRW) models whose cosmological fluid 
possesses a spatially homogeneous and isotropic refraction index n(t). We use the optical 
metric to derive both apparent-size, dA(z n ), and luminosity, dL(z n ), distance-redshift rela- 
tions. In Section IVT1 we fit the di(z n ) of our index of refraction model to current supernova 
data. Without a cosmological constant A, we also estimate the value of n(t) required to 
produce roughly the apparent size distance d^iz) at last scattering required by the position 
of the acoustic peaks in the angular power spectrum of the WMAP data. In Section l"VHI we 
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draw our conclusions. 



II. THE OPTICAL METRIC 

In GR type theories, a gravity field is described by a metric g ab on a four dimensional 
manifold (we use a +2 signature here). We additionally assume the presence of an arbitrarily 
moving medium with normalized 4- velocity u a u a = — 1 that fills spacetime. For simplicity, 
we also assume that the fluid's electromagnetic properties are linear, isotropic, transparent, 
and non-dispersive; and can be summarized by two scalar functions: a permittivity e(x a ) and 
a permeability fi(x a ). Following Ehlers [3| in this section, we write the two electromagnetic 
bivectors as F ab (B,E) and H ab (H, D). They satisfy Maxwell's Equations 1 

d[ a F bc] = 0, 

= — J a , (1) 

c 

with constitutive relations 

H ab u b = eF ab u b} 

F[abU c \ = fiH [ab u c] . (2) 

(For a familiar example, take Minkowski spacetime with the fluid at rest, u a = (1, 0, 0, 0).) 
The optical metric of Gordon [1] is defined as 

g a b = g a b + (l - — )u a u b , (3) 
efi 

with inverse 

g ab = g ab + (1 - ^)u a u b . (4) 

The combination efi is related to the usually defined refraction index n(x) = y/eJI. To relate 
covariant derivatives of the two metrics and to obtain Eq. ([9]) below the relationship of the 
two determinants is needed 

det{g ab } = —det{g ab }. (5) 
efi 

At this point we have one differentiable manifold with two metrics or equivalently two 
related spacetimes, the physical and the optical. In this paper we are primarily interested in 



1 Square brackets [ ] symbolize complete anti-symmetrization of the enclosed indices. 
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the dynamics of a particular type of physical field (radiation) in the optical spacetime. All 
physical objects are described by tensor fields in physical spacetime and those of interest here 
have associated fields in optical spacetime. Where necessary we denote optical spacetime 
fields with a bar. The optical equivalent of the physical covariant Maxwell field F a b is F a b 
itself; hence the homogeneous Maxwell equations are satisfied in both spacetimes, and both 
share covariant 4-potentials. Using the optical metric the two constitutive equations can be 
written as a single equation 

H ab = -F ab = -g ac g bd F cd , (6) 
/i /i 

as can be seen by first expressing the contravariant Maxwell field in optical spacetime as 

F ab = [F ab - (1 - efi)u a F b c u c + (1 - tii)u b F a c u% (7) 

and then contracting separately with Ub and e a b c dU c , where e a b c d is the completely anti- 
symmetric Levi-Civita symbol. Here the metric dependent bivector F ab is the metric inde- 
pendent 2-form F a b raised using the optical metric Eq. (J4]) rather than the physical metric g a b- 
By using the identity relating the contracted Christoffel symbols to the metric's determinant 

{^} = <9 c lo gv /-det{g ab } , (8) 

and Eq. ([5]), Maxwell's equations ([1]) can be written using the optical metric as 

d[aF bc ] = 0, 

V fe (e 2 *F ba ) = —J a = —JqiJ\ (9) 

v ' c c 



where e 2(/i = y e//i and the covariant derivatives are taken using the optical metric (see Ehlers 
[if] ] for details when J a = 0). The form of the inhomogeneous equation is slightly modified 
and hence slightly more complicated, i.e., the e 2 ^ term is present in Eq. (Q; however, the 
advantage is that there are no constitutive equations ([2]) to deal with, i.e.,F ab is just F ab 
raised with the optical metric. Solutions constructed in the optical spacetime via Eq. (Q 
directly translate to solutions in physical spacetime via Eq. (jH]). Because the vacuum Maxwell 
equations are conformally invariant, any metric conformally related to Gordon's can be used 



to generate H ab , see 16[. In the next section we show that light waves travel along null 
geodesies at the speed c in the optical spacetime (and hence in any conformally related 
spacetime), whereas the corresponding waves travel at speed c/n in physical spacetime. 



III. GEOMETRICAL OPTICS APPROXIMATION 



In this section we follow Sachs [15( and Ehlers [17| but use vectors rather than spinors 
to derive the transport equations for the amplitude and polarization of a geometrical optics 
wave. Readers not interested in the tensor calculus details can skip this section. We assume 
that the electromagnetic wave is planar on a scale large compared with the wavelength, but 
small compared with the curvature radius of spacetime. We write the covariant (and metric 
independent) field tensor as 

F ab = H |e 45 / A (a* + *B ab + 0(A 2 )) | , (10) 

where A is a wavelength related parameter, S(x a ) is the so-called eikonal function and is 
real, and 3?{-} stands for the real part. The A ab term represents the geometrical optics 
(GO) approximation and the B ab term is its first order correction in both the physical and 
optical spacetimes. Defining the unitless (also metric independent) wave vector k a = d a S 
and inserting Eq. (11 01) into the vacuum Maxwell equations (J a = in Eq. (Q) we obtain to 
order A -1 

A[abk c ] = 0, 

A ab k b = 0, (11) 

and to order X° 

d[ a A bc } + k[ a B bc ] = 0, 
V b A ab + B ab k b + 2A ab (f) >b = 0. (12) 

All barred contravariant quantities throughout are obtained by raising indices with the 
optical metric, e.g.,A ab = g ac g hd A C( i\ unbarred are obtained by raising with the physical 
metric g ab . Equations ffTTj) tell us that k a = g ab k b is tangent to null geodesies of the optical 
metric 

k b V b k a = 0, (13) 



and that A ab is of the form: 



A ab = -2k [a £ b] , (14) 
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where £ a is spacelike and constrained by £ a k a = but has the remaining freedom of definition 
£a —> &a + f{x)k a . It is £ a that determines the amplitude and polarization of the GO wave 
seen by an observer and it is Eqs. (TT5|) that establishes the speed of propagation as c. 
Equation f[T2"j) tells us that 

B ab = 2(£[ a , 6 ] - k[ a T>b]), (15) 

with a remaining freedom V a — > V a + g(x)k a and gives as the propagation equation for £ a 

+ + £ a = y (V 6 £ 6 + A; b P b + 20, 6 £ 6 ). (16) 

The affine parameter derivative symbolized by ' ' ' is the invariant derivative Ar V& along the 
null geodesies generated by the vector field k a . If we now split £ a into a scalar amplitude 
and a unit polarization vector, i.e., 

£ a = £e a , £>0, e a e* a = 1, (17) 

where * means complex conjugate, the transport equation for the amplitude S becomes 

£ + £9 + 8<j> = 0, (18) 

where 9 is the expansion rate of the null rays defined by the vector field k a . It is defined by 
the divergence of k a 



9 = -V a k a = ^L, (19) 
2 JA K J 



and is related to the fractional rate of change of the observer independent area A of a small 
beam of neighboring rays [l^. Given A, we are able to integrate Eq. (1181) 



£ Ef^tef^' <2o) 

where the subscript e means evaluate at (or close to) the emitter. 

For the calculation at hand we need the amplitude £ only, however, if we were interested 
in the wave's polarization a suitable choice for f(x) makes the right hand side of Eq. f|T6|) 
vanish and also makes e a = 0, i.e., a particular choice for a polarization vector can be made 
that is parallelly transported along the null geodesies of the GO wave. 

The GO approximation is just the 0(X°) term in Eq. (fTOl) i.e., 

F ab = -2£ ^{e iS/x k [a e b] }. (21) 



The frequency and wavelength seen by observers comoving with the fluid can be computed 
using the fact that the phase of the wave changes by 2ir when the observer ages by one 
period of the wave r, or respectively steps a spatial distance of one wavelength A in the 
direction of the wave, k a , 



2tt 



S 



A 



CT 



CTU a )d a (-)=--(u a k a 



27T = (Xk a )d a [j 



A 

~n (u a k a ). 



(22) 



The natural choice of the constant parameter A is the rationalized wavelength (A e /27r) at 
the emitter. This requires the eikonal satisfy n e (k a u a ) e = —1. 

The energy and momentum of this wave as seen by a comoving observer in physical 
spacetime {u a g ab u b = —1) is contained in the Poynting 4-vector 



S" 



C 



H ac p 



cb 



\s a b H dc F cd 



u 



(23) 



where all quantities in Eq. (J231) are evaluated in the physical spacetime. When evaluated 
using Eqs. (El) and ( 12TT) 

1 



S a 



F ac F 



cb 



:S a b F dc F cd 



u 



S 2 k a (k b u b ). 



(24) 



In the last equality the oscillations have been averaged over. The energy density and 3-d 
Poynting vector seen by observer u a are respectively 



U = --S a u a 

c 



-E\k a u a ){k b u h 



n 



57T/i 



-S 2 (k a u 



a\2 



s a 



S a -cUu a 



inn 



-£ 2 {k b u b ) k a + n 2 {k c u c )u a 



with magnitude 



S i = \ S a \S 



en 



>_La 



inn 



8 2 {k a u 



a\2 



A P T% 
A T 2 ' 



C - 1 e 'e 

<->± e — r — ■ 



(25) 



(26) 



In the last equality we have eliminated the amplitude £ using Eq. (|20l) and continue using 
a subscript e to represent quantities evaluated near the emitter. Equation ( |26i) simply says 
that the energy flux varies inversely with the beam's area and inversely with the square of 
the period, even in the presence of an index of refraction. We will use this expression in 
the next section to compute the luminosity distance-redshift relation for FLRW cosmologies 
that are filled with a transparent optical material. 



7 



IV. OPTICAL METRIC FOR ROBERTSON- WALKER (RW) SPACETIMES 



The familiar spatially homogeneous and isotropic Robert son- Walker (RW) metric can be 
written as 



ds 2 = _ c 2 dt 2 + R 2,fi [ + r 2< dQ 2 + gin 2 q 1 (2?) 

[ 1 — «H J 

where k = 1,0,-1 for a closed, flat or open universe, respectively. The cosmic fluid as- 
sociated with the RW metric is at rest in the co-moving spatial coordinates (r, 8, 0) and 
hence has a 4- velocity u a = 5f/c. We assume that the cosmic fluid has associated with it a 
homogeneous and isotropic refraction index which depends only on the cosmological time t, 
i.e., y/ejl = n{t). From Eq. ([3]) only the g tt component of the optical metric is seen to differ 
from the physical metric, i.e., 

ft< = ■ (28) 

The radial null geodesies of the optical metric are found by fixing (9, 0) and integrating 

Two such geodesies traveling between the origin and a fixed comoving point r satisfy 
H ° cdt rt +At C( fc rr dr 



lt e n(t)R(t) Jt e +At e n(t)R(t) Jo y/l - kr 
where we have defined 

{sin[r] k = +1, 
r k = 0, (31) 

sinh[r] k = — 1. 

In the above equation, (t e ,t Q ) and (t e + At e ,t Q + At D ) represent the emitting and receiving 
times of the two respective null signals. We see from Eq. (130]) that the differences in emission 
and observation times are related by 

At, At, 



el 



(32) 



n(t )R(t ) n(t e )R(t 
and hence the redshift z n is given by 

l + z n = ^= n{ ' l R{ '°l =^(l + z). (33) 
At e n(t e )R(t e ) n(t e ) y ' 1 ; 

We have used z n as a measure of the observed frequency change but have also kept the usual 
z which measures the wavelength change. The tangent to the radial null geodesic (k a k a = 0) 
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can be found directly from Eq. (|29|) by a variation 



n y/1 — kr 2 



k "= a Ufl'^^' ' ■ (34) 



with covariant components 



ka = a \^R' T/T^W' °' °) ' (35) 

The constant a is arbitrary and equivalent to the freedom of choosing an affine parameter 
along a null geodesies. In the physical metric the light ray has a timelike tangent vector, 
i.e., 

k a g ab k b = n 2 (l - n 2 )(k a u a ) 2 = a 2 ) < 0. (36) 

The eikonal S of the GO approximation Eq. (1211) can easily be found for this covariant 
vector field k a assuming the spherical wave originates from an emitter located at r = 

To relate a and the constant A of Eq. (|2ip to comoving wave length, we use Eqs. (1221) and 

The last equality results from choosing a = R(t e ) and A = X(t e )/2ir and confirms our 
interpretation of the conventional (1 + z) — R/R e as the wavelength redshift, see Eq. fl22l) . 

We are interested in computing the apparent-size and luminosity distances for RW models 
with an index of refraction n(t) and must be careful in doing so. The optical metric gives 
the correct wave trajectories, but because it does not measure distances or times correctly, 
densities and rates will be incorrect. Because angles, areas, and redshifts are easier to 
calculate than energy fluxes we start with the apparent size distance-redshift d,A_{z n ). We 
will then compute the luminosity distance di(z n ) by using the 3-d Poynting vector of Eq. ( |25i) . 
We use the optical metric in the form given in Eq. (1281) with Eq. (12 7p because the coordinates 
(t, r, 6, 4>) have direct physical interpretations in the RW metric itself. For an example, in 
the local rest frame of the source (observer), the proper time interval is At e (respectively 
At a ) instead of At e /n(t e ) (respectively At /n(t )), and hence the observed shift in periods, 
At Q / At e , is correctly given by Eq. (|33|) . From Eq. (1271) the apparent size distance (also called 
the angular size distance) of a source at coordinates (t e , r) is just 

d A = rR(t e ), (39) 



as seen by an observer at (0, t Q ) where the three coordinates (r,t e ,t ) are constrained by 
Eg. fl30|) . To give dA(z n ) we must use Eqs. ( 130]) and fl33|) to eliminate (r, t e ) in terms of 
(z n , t Q ). We start by using Eq. (I3"3l to change variables in the remaining integral of Eq. (I3"U|) 
from t to z n . The following steps are familiar except for the presence of the index of refraction 
n(x) and the two redshift variables (z n , z). The dynamical equations of Einstein are used to 
change from t to R and then to (1 + z) = R Q /R. 



sinn 1 (r) 



*° cdt f R v cdR 



, n(t)R(t) Jr £ n(R)R(dR/dt) 

c r z ( z ^) dz 



(40) 



where 



and where 



RqHq Jo n(z)h(z) 

h(z) = y/[n A + n k (i + zf + n m (i + ?Y + n r (i + z y], (4i) 



n k ee --^-2 = i - (o A + a m + n r ). (42) 



c 2 fc 

n n 

The wavelength redshift z(z n ) as a function of the frequency redshift z n is found by elim- 
inating t e in Eq. (1331) . We refer to solutions to Einstein's equations with a RW symme- 
try as Friedmann-Lemaitre-Robertson- Walker or simply FLRW cosmologies. The three Q 
constants represent, as usual, current relative amounts of non-interacting gravity sources: 
vacuum, pressureless matter, and radiation energies. From Eq. (I3"$j) we conclude that the 
apparent size distance-redshift relation for a FLRW cosmology with an index of refraction 
is 

, , , 1 c 1 " r— 7 r z ( z ") dz 

d ^ = (TT7R))^^v sum '■ ' (43) 



/o n(z)h(z) 

To derive the luminosity-redshift relation diizn) knowing dA(z n ) one ordinarily uses 



Etherington's 18| result 

d L (z) = (l + z) 2 d A (z). (44) 

If this result were correct with an index of refraction present, one would need to know 
which redshift to use, frequency z n or wavelength z. To know what to choose we evaluate 
the magnitude of the Poynting vector Eq. (1261) and arrive at the correct replacement for 
Eq. (J44]). To find the needed area A we evaluate the expansion 6 of Eq. (fT9l) using Eqs. (|34j) . 
(J27D, (TJHD, and ©. We find a simple result 
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and hence the beam area A oc rR from which we have the needed Poynting vector magnitude, 
Eq.flM]), 

4vr(r^ r e 2 _ L e _ L e 

1 ±e 4<n{rR) 2 T 2 4<Il{rR) 2 {l + Z n ) 2 47id 2 L y ' 

The latter identity defines the luminosity distance di in terms of the total power radiated 
at the emitter L e and the flux received, i.e., the Poynting vector at the observer, 

d L = rR{l + z n )=rR e (l + z)(l + z n ), (47) 

which agrees with the Etherington result Eq. (1541) only if we use one frequency redshift factor 
(1 + z n ) from Eq. fl33|) and one wavelength redshift (1 + z). Equation (j4"6]) also confirms a 
conserved photon number interpretation of the radiation even in the presence of a time 
dependent index of refraction (which has the potential of taking energy out of the radiation 
field) . It is equivalent to having a fixed number of photons emitted in a time At e each having 
energy hv e and all being collected over an area 4n(rR) 2 in a time At Q but with redshifted 
energy hv Q . 

V. AN EFFECTIVE INDEX OF REFRACTION INDUCED BY THE SACHS- 
WOLFE EFFECT 

Up to now, we have been considering a refractive index modeled after the one generated 
by induced electromagnetic polarizations in inter and intra galactic media, dark or otherwise. 
Lensing has long been interpreted as a gravitationally induced refraction effect, and here 
we suggest that to I s * order, inhomogeneities in the flat FLRW models are equivalent to 
effective indices of refraction. 

Sachs & Wolfe 



25j ] were two of the first to consider the effect of perturbations of the 
homogeneous and isotropic models on optical observations. In that classic paper, the authors 
used perturbations in the flat, i.e., k = 0, FLRW spacetime to study the angular fluctuations 
in the CMB. They used a conformally flat version of the metric 

ds 2 = R 2 (r]) [r] ab + h ab ] dx a dx b (48) 

with dimensionless coordinates and worked in a comoving gauge to derive the equations 
governing the evolution of the metric perturbation h ab and perturbations of the energy- 
momentum tensor 5T ab . Here the conformal time coordinate of the flat Minkowski metric 
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rjat is x° = 7] and for the pressureless case is familiarly related to the comoving FLRW 
time coordinate t by 77 = (3 tHo/2) 1 ^ 3 . The three Euclidean spatial coordinates are labeled 
by letters of the Greek alphabet. They then considered the deviations of null geodesies 
from the unperturbed case and derived the now famous temperature fluctuations in the 



micro-wave background caused by h ab} see Eq. (42) of Ref. [25j. Among the scalar, vector, 
and tensor perturbation modes in dust (p = 5p = 0, see Eq. (22) of Ref. [251]) . the scalar 
density perturbations, i.e., the relatively decreasing A(x J ) mode and relatively increasing 
£>(x 7 ) mode (responsible for the famous Sachs- Wolfe effect) give 

1 rf 

h a p = ——^A ;a i3 + 5 a pB + —B !(X p , (49) 

and h 0a = 0. The arbitrarily specified form of the scalar modes are related to the density 
perturbation 5p through Poisson's equation 

H% 9 f6A 3B\ 



(See Eq. (22) of Ref. 25j). The time component of the perturbation h 0Q vanishes because of 
a comoving gauge choice (u a oc 8q) and ho a is only present for the rotational perturbations. 
To connect this metric to Gordon's optical metric we must make the following non-gauge 
change of coordinates 

x a = x a + —^A !a (x) — —B >a {x), 
V = v(l- ^A(x) - ±B(x^j , (51) 

and obtain to 1 st order 



ds 2 = IIH'I) (l " ^A(x) + ^B(x)) \ ^ : dr ■:■ r-UW- : sm- 

\ t 10 J \ l-S A (x) + lB(x)) 



(52) 

This simply says that the Sachs- Wolfe metric is a conformally transformed Gordon met- 
ric corresponding to a k=0, FLRW metric with a spacetime index of refraction n = 
1 — 6A(x)/f] 5 + 3B(x)/5. Since conformal transformations don't alter light cones (see [if]]) 
we have arrived at the connection of null geodesies of Sachs- Wolfe's density perturbations 
with the light curves of an unperturbed FLRW spacetime possessing an index of refraction. 
In contrast to the homogeneous optical fluid discussed in Section HVl the comoving frame 
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of the index of refraction in Eq. (1521) is not the same as the comoving frame of the matter 
density in Eq. fj48|) . However, they are related by the coordinate change of Eq. (15"TT) . 

We have our doubts about extending the index of refraction comparison beyond linear 
perturbations, and make no claims as to that possibility. Such an extension would be quite 

on non-linear observational effects in Swiss Cheese 
IE 2! lj, also hoping to find sources of apparent 



interesting because old work [26 
cosmologies are again in the literature 



acceleration other than a cosmological constant. Work on interpreting effects of local density 



perturbations on the Hubble curve are numerous 
of which can be compared to the above in the 1 st order regime. 



mm 

, [11I.I12I. [13I . Q, the results 



VI. FITTING SUPERNOVA DATA WITH A REFRACTION INDEX MODEL 

In this section we use the index of refraction model of Section [IV] to fit the current 

32]]. We use the 178 supernova from the gold sample 2 with 



supernova data [28 



29, 



30 



31 



redshifts greater than cz = 7000 km/s , see Fig.HJ The Hubble constant we use is Hq = 65 
km/s/Mpc and since we are concerned with the matter dominated era, we exclude radiation 
(Q r = 0). We compare the distance modulus versus redshift, n(z), of the concordance model, 
— 0.7, Q m = 0.3, n = 1 with two n(z) > 1 models. The first is a baryonic matter only 
model (Q m = 0.04) and no cosmological constant (Q\ = 0) with n(z) = 1 + O.lz 2 — 0.0452; 3 . 
The second model includes a dark matter contribution, Q m = 0.3, no cosmological constant, 
and n(z) = 1 + 0.15z 2 — 0.05z 3 . Also included is a now disfavored dark matter only model, 
= 0, Q m = 0.3, n — 1. In the inset of Fig.[T], we use this case to compare with the 
two n 7^ 1 models and with the concordance model. The critical redshift region is between 
0.2 < z < 1.2, where most of the supernova data is concentrated. Both n 7^ 1 models 
fit the data much better in this region than models with the same Q parameters but with 
no refraction. The two refraction indices are plotted in the insets of Fig. [21 As the reader 
can easily see the effects of a suitable index of refraction n(z) can simulate the accelerating 
effects of a cosmological constant. 

The Supernova data is currently considered the best evidence for the existence of dark 
energy because of the observed acceleration in the expansion of the universe (see Fig. 7 of 



http://braeburn.pha.jhu.edu/~ariess/R06/. 



13 



Ref. (29||). For the homogeneous FLRW models, the acceleration is directly related to density 
and pressure by 

(53) 



R 
R 



-(p + 3-) 



3 - 

A true observed acceleration, i.e., R > 0, requires p < — pc 2 /3, and hence implies an unusual 
equation of state such as vacuum energy (p = —pc 2 ). What we show here is that an over- 
looked index of refraction can cause a misinterpretation of the Hubble curve, suggesting an 
acceleration. In Fig.[2]we plot H(z) and R(z)/R Q = H(z)/(1 + z) for the two n^l models. 



They can be c omp ared with similar plots in Ref. 



flux averaging 33 



291 ] . The data points are plotted using 



341 ] and uncorrelated redshift binning 35] algorithms. To apply these 



techniques to non-flat cases, we define 

rz(Zn) 



g{Zn) = 



dz 



o n(z)h(z) 



=smn 



1 + Z n 



101 



(54) 



where \i is the distance modulus 



li — 5 log 



1 Mpc 



25. 



We furthermore defined 



<?(4 +1 ) 



(55) 



(56) 



£ i+1 (4 +1 )-^(4)' 

which when averaged inside each bin gives us an estimate of the inverse of the product of 



n(z) with h(z) = H(z)/Hq (compare with Eq. (5) of Ref. (35j). The presence of an index of 
refraction produces a degeneracy in determining the value of H(z) and hence in R(t). A 
suitably decreasing n(t) and a non-accelerating R(t) will mimic an accelerating universe. 

We can see that our index of refraction models fit the data well. However, we need to 
remind the reader that the binned data plotted on the H(z) and R(z)/Rq curves are model 
dependent. The binning process as designed in Ref. [35(] requires use of di(x n ), i.e.,g(z n ) in 
Eq. (I54p . Rather than using this technique to argue for an observed accelerating H(z), we 
argue for an observed n(z) with a non-accelerating H(z). The point we make is that we can 
fit the /i(z) data with no A, and are able to get rid of the acceleration. 
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FIG. 1: Distance modulus /i versus redshift z. Dashed curve: Oa = 0,$7 m = 0.3, n = 1; Green 
curve: Q A = 0.7, Q m = 0.3, n = 1; Red curve: ft A = 0, fi m = 0.04, n(z) = 1 + O.lz 2 - 0.045z 3 ; 
Blue curve: Qa = 0, Q m = 0.3, n(z) = 1 + 0.15z 2 — 0.05z 3 . A^(z) are given in the inset. The 
difference is taken with respect to the fiducial case where f^A = 0, O m = 0.3, n = 1. 




FIG. 2: H(z) (upper panel) and R{z)/R Q (lower panel) curves for the two n / 1 models. Left 
column parameters: $7 a = 0, f2 m = 0.04, n(z) = 1 + O.lz 2 — 0.045z 3 ; Right column parameters: 
n A = 0,n m = 0.3, n{z) = 1 + 0.15z 2 - 0.05z 3 . 

VII. FLATNESS OF THE UNIVERSE 



The conclusion drawn from the latest WMAP data, when combined with the SNe la Hub- 
ble curve, that vacuum energy exists, depends crucially on many unconfirmed theoretical 
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500 



1000 

z 



1500 



FIG. 3: Dashed curve: f^A = 0, J7 m = 0.3, n(z) = 1; Blue Curve: Oa = 0, J7 m = 0.3, n(z) 
1 + 4.45 x 10 - V - 1.25 x 10" 10 z 3 ; Green curve: ft A = 0.7, VL m = 0.3, n = 1. 



assumptions including the adiabatic power law assumption for the initial perturbation spec- 



trum [36J, 137J]. Such observations have motivated efforts to find ways to produce a perceived 
acceleration other than by a real A , |l0j, lll|, ll2|, ll3|, Il4| . This section is another 

such effort. 

The angular position of the first acoustic peak is commonly believed to be the strongest 
piece of evidence for the flatness of the universe. The characteristic wavelengths of the 
acoustic oscillations at the last scattering surface depend very weakly on A, but their ob- 
served angular size as seen by us now depends significantly on a combination of Q m and ^a, 
see Eq. (143 p . Assuming our Universe is of the FLRW type with no refraction index, a first 
acoustic peak at ~ 0.8° is almost fit by a flat universe. 3 With a suitable index of refraction 
and no cosmological constant we can produce an angular diameter distance comparable to 
the angular diameter distance of a flat cosmology at any given redshift, independent of the 
Hubble parameter H . In Fig. [3] we have used a Q\ = 0, Q m = 0.3 model with an index of 



3 Within the context of a power law ACDM model (w = 1), WMAP data alone does not rule out non-flat 
models. With a prior on the Hubble constant, Ho > 40kms -1 Mpc~ 1 , or combined with other astronomical 
observations, such as SDSS LRG sample, HST constraint on the Hubble constant, or SNe data, WMAP 
data strongly favors a nearly flat universe with nonzero vacuum energy, see Table 12, Fig. 20, and Fig. 21 
of 37|. For a more general model of dark energy, e.g., one with a time evolving equation of state parameter 



m^l instead of a cosmological constant A, significant spatial curvature is still allowed even when Hq is 
not restricted to be small, see e.g. [38 
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refraction n(z) shown in the inset. We chose this n(z) because it produces similar distances 
over the large redshift range z > 1000. To conclude that WMAP implies flatness requires 
the acceptance of the accuracy of theoretical assumptions beyond the initial perturbation 
spectrum; e.g., even the accuracy of the optics of homogeneous FLRW models is now being 
questioned as was pointed out in Section [V] 



VIII. DISCUSSION 

We have reviewed Gordon's optical metric theory [l| which incorporates an index of 
refraction into its geometry. We then used the optical spacetime to derive the transport 
equations for the amplitude and polarization of a geometrical optics wave. We applied it to 
the homogeneous, A = 0, FLRW models and estimated the refraction index needed to fit 
current SNe la and WMAP data. We found that an n(z) ~ 1.07 at redshift z = 1.5 in a 
baryon only model, or an n(z) ~ 1.15 at z = 1.5 in a dark matter model, could easily fit the 
supernova data (see Fig.fU and Fig. [2]), and that an n(z) as big as 1.3 at the last scattering 
surface in a dark matter model would give the same angular diameter distance (Ia(z) as the 
concordance model (see Fig. [3]). 

The question is, where could such an index of refraction come from? If it had its origin in 
atomic dipole moments or charges in plasmas the densities would have to be much larger than 
they actually are. A critical mass density now is about 8 x 10~ 30 g • cm -3 , which translates 
to 1 x 10~ 20 g-cm -3 at z = 1100. The density of air on the earth is about 1.2 x 10~ 3 g-cm~ 3 , 
some 10 17 times denser than the universe at recombination and yet its index of refraction is 
only n = 1.0003. We conclude that there is little hope for a baryon- lepton origin for n. A 
long shot would be a colorless index of refraction for the mysterious dark matter. 

Severe limits have already been estimated on direct interactions of the dark matter 



particles with photons caused by fractional charge 40j, |4l| (q/e < 10 -5 — 10~ 7 depend- 
ing on the particle's mass) and by electric/magnetic dipole moments 42| (dipole moment 
< 3 x 10~ 16 ecm). General limits on photon interaction rates have even been estimated by 
requiring the associated collisional damping scale be small enough to allow structures larger 
than 100 kpc to form 43(. Proposing an n of unknown source is perhaps outlandish, but 
not much more than proposing a non-intuitive repulsive cosmological constant A to produce 
acceleration. Even though the latter has become fashionable, we wish to add a dark index 
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of refraction theory to the lists of alternatives to think about. 

In this paper we developed the general framework needed for using Gordon's optical 
metric in cosmological observations, but have applied it only to an index of refraction model 
which is homogeneous and isotropic. If the dark matter and it's assumed index of refraction 
were truly homogeneous we could have additionally proposed a redshift dependence for n(z) 
modeled after a dilute dielectric gas or plasma. However, such a model would still contain 
unknowns equivalent to ionization densities and/or molecular polarizabilities. Instead we 
chose a phenomenological expression in the form of a cubic containing two parameters which 
we adjusted (i.e., n(z) = 1 + az 2 + f3z z ). Such a simple starting point is prudent because 
we know the real universe is filled with low density voids, and high density condensations, 
as well as associated velocity perturbations all of which would modify the refraction index 
n{x a ). If an index of refraction model such as the one proposed here has merit, future 
efforts can look into how such perturbations, including local variations in the magnetic field 
of the intergalactic medium, might impact distance- redshift. However, the optical metric 
theory Eq. (jHJ) is still the applicable theory. We also leave to the future, further exploration 
of the equivalence of the optical effects of gravitational inhomogeneities (beyond the linear 
perturbation results of Sachs- Wolfe in Section |V|) and our index of refraction proposal. 
Complete equivalence would be quite interesting and useful in light of the current interest in 



Swiss Cheese optics |2|, I3J, 0, 2y, |27j. Modifications in distance-redshift caused by random 



spacetime perturbations could then be interpreted as being caused by an effective index of 
refraction. The idea of an optical metric can also be applied to other massless particles 
which follow null geodesies in vacuum. If the presence of material causes interference of 
the propagating waves of the particles, an effective refraction index should exist. For an 
example, gravitons and some neutrinos are massless and local inhomogeneities, such as the 
Sachs- Wolfe perturbations discussed in Section [V] would alter propagation of their waves 

The idea of solvin g co smological problems via a changing light speed model is not new, 



sec 



44 



45 



4f|. Ellis 47] has recently pointed out consistency constraints required of such 
theories. Our proposal is fundamentally different from those cited above in that it is based 
on a classical electrodynamics analogy (the cosmological fluid simply has an unexpected 
refraction index which reduces the propagation speed of electromagnetic waves). Because 
we are not proposing a change in the vacuum light speed c or the limiting causal speed, the 
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proposal is not subject to Ellis's criticism. 

Finally we note that since an accelerating universe is consistent with other observations, 
such as Baryon Acoustic Peaks detected in galaxy surveys 



49 



50 



5l| and the interesting 



H(z) relation obtained from ages of passively evolving galaxies in [52j|, additional compar- 
isons with refraction models are in order. 
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